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Abstract 

The underlying mechanism for Adaptive Feedback Control in the experimental photoisomerization of 
NK88 in methanol is exposed theoretically. With given laboratory limitations on laser output, the compli- 
cated electric fields are shown to achieve their targets in qualitatively simple ways. Further, control over the 
cis population without laser limitations reveals an incoherent pump-dump scenario as the optimal isomer- 
ization strategy. In neither case are there substantial contributions from quantum multiple-path interference 
or from nuclear wavepacket coherence. Environmentally induced decoherence is shown to justify the use 
of a simplified theoretical model. 

PACS numbers: 82.50.Nd, 82.53.Uv,82.30.Qt 



1 



Experimental results on the quantum control of molecular processes yj, |2J] fall into two cate- 



gories: those designed to exp 



ore the utility of a particular coherent control scenario, and those that 



use adaptive feedback to attempt automated optimal control of a target process. 

The former have well understood mechanisms but are scenario-specific, whereas the latter are 
generally applicable but have thus far provided very limited insight into the mechanisms by which 
control is achieved. Indeed, it is widely recognized that extracting the mechanism of control from 
an optimal control experiment is the central challenge in this research area [9]. In this letter we (a) 
expose the simple mechanism underlying control in a recent [7] trans-cis isomerization control 
experiment, and (b) demonstrate that removing experimental restrictions on laser frequency and 
amplitude exposes modified versions of well known coherent control scenarios as the dominant 
control mechanism. 

Gerber's seminal experiment |Q] shows successful optimization of the trans to cis isomerization 
of 3,3'-diethyl-2,2'-thiacyanine iodide (NK88) as well as suppression of the trans to cis transition, 
depending upon laser pulse shapes. As is typical of these experiments, the optimal pulse shapes 
in both cases are complex functions of phase, frequency and time, and differ substantially for the 
two targets. Below we obtain quantitative agreement with these results. Further, we do so, despite 
the complexity of the system, using a one degree of freedom model coupled to a bath. Indeed, 
it is precisely the presence of environmentally induced decoherence that simplifies the control 
problem. 

Consider a model consisting of the system with Hamiltonian Hs, the bath H B , system-bath 
coupling Hsb, and system-electric field coupling described in the dipole approximation. The total 
Hamiltonian is given by 

H = H s +H B +H SB -LiE(t), (1) 

where Hs describes the isomerization process via a one dimensional reaction coordinate <j), H B 
represents all other degrees of freedom, /i is the transition dipole moment, and E(t) is the incident 
electric field at time t. 

In terms of the two participating electronic states, the system Hamiltonian is given by: 

\ Ve g (<f>) K + V^)J 

where K = — j^jp is the kinetic energy, Vg(0) and V e (0) are the ground and excited electronic 
state potential surfaces, and V r ge (^)) = V e g(0) is the coupling potential between ground and excited 



states. In the adiabatic representation the ground state potential is a double well fl. Populations 
reported below are of the eigenstates of the full molecular Hamiltonian H^. The simplest dynamics 
takes place by photoexcitation from the trans configuration to the excited electronic state followed 
by de-excitation to the cis and trans ground state via system-bath coupling. 

In the case of NK88, the effective mass is m = 5 amu A 2 , V g (0) = Aq (1 — cos0), V e ((j)) = A\ + 
A 2 cos0, V eg (» = 1000 cm" 1 , where Aq,A\ andA 2 are 15900cm- 1 , 17500cm- 1 and 7500cm- 1 , 
respectively. These, and other parameters below were obtained by a fit to experimental results. 
The eigenvalues A,- and eigenvectors of the system are calculated by diagonalizing the molecular 
Hamiltonian matrix represented on a grid, where periodic boundary conditions at 2n are imposed 
on the system. The transition dipole moment ji is expressed in terms of the two electronic states 
as a 2 x 2 dimensional matrix with zeroes on the diagonal, and off-diagonal elements jU ge (0) = 
lODebye. This corresponds to the oscillator strength / ~ 1 in the Franck-Condon region of the 
trans configuration. 

The bath is described as a set of harmonic oscillators of frequency co a and the system-bath 
coupling is H^b = QUa nK a (ba + b a j , where b\ and b a are the creation and annihilation opera- 
tors pertaining to the ath harmonic oscillator. The operator Q is a diagonal 2x2 matrix with cos <j) 
on the diagonal, and the coupling constant K a and spectrum of the bath are chosen in accord with 
an Ohmic spectral density /(<») = 27iX a K 2 <5((y — (O a ) = r\(De~ a l ac , where the strength of the 
system-bath coupling is determined by the dimensionless parameter r\ = 5, and co c = 450 cm- 1 . 
Given these parameters, the electronic dephasing time around the Franck-Condon region of the 
trans configuration is ~ 10 fs, and virtually complete relaxation from excited trans to stable trans 
and cis occurs within 5 ps. The former is a typical characteristic dephasing time whereas the latter 
is chosen to agree with experiment. 

The dissipative dynamics of the system was evaluated using the Redfield equation with secular 
approximation lioll . The evolution of diagonal elements Pu(t) of the system density matrix is 
given by 

i-Pii (t) = -V^p- V [pimit) jlmi ~ fMmpmi{t)] 

+Y,wijpjj{t)-pu(t)Y,wjh (3) 

where the transition probability is Wji = Ffj .. + T^.. and where each index denotes a state of the 
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system, including the electronic and vibrational quantum numbers. Here, T^ ik — i ± 



r+., and 



1 



Hijl = ^QljQik L dT I d(OJ(CQ 



2n - jo Jo 
[n(co) + 1] e -'^ ik+C0)T + n(co) e'^-^ 



(4) 



where n(co) = {exp(ftco / k^T) — 1} 1 is the Bose distribution, T is a temperature, and COji = 
(Xj — Xn /h. After some algebra, we obtain 



Wji 



2 

\Qji\ J (-(Oji) [n(-(Oji) + l] for cojj < 



\Qj i fj((Oj l )n((O jl ) 



(5) 



for COji > 



The evolution of the off-diagonal elements is described as 



— Pij ( t ) = -i(OijPij{t) - Yijpij(t) 
E(t 



i j_ \ Pimjt^) fim j H-impmjit 



with dephasing rate yn: 



jkkj 



r r . 
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The resultant vibrational dephasing time within the excited electronic state is « 15 fs. 

To model the adaptive feedback experiment, the electric field comprises 128 frequency values, 
where the phases of each frequency component are the optimization parameters; the frequency 
width is 200 cm -1 , and the time width is 2ps. Specifically, the electric field function is therefore 
taken to be 



E(t) 



127 

E^exp 

;=0 



• cos (£ty + 0,-) 



to 



2At 



Zlj — Q.q 
2A£l 



(8) 



where A=5 MV/m, t = At = 2ps, C2 = 25000 cm" 1 , = 200 cm" 1 , Q, = 24800 + 3.125 x i in 



cm x , and 0/ are the optimization parameters. Note that in accord with experiment [7] the field 
is optimized by varying the phases ; using an evolutionary algorithm, and the field amplitude, as 
well as the overall frequency width of the pulse, are constrained. Further, the algorithm is designed 
to simulate experimental conditions I 111 , where the population size is 60, 10 survivors are selected 
from the generation, each of which has 4 children by mutation and 1 child by crossover, and 
yielding again 60 individuals for the next generation. 



In our second study, the experimental frequency and amplitude constraints were lifted and 
optimal field was evaluated by using a monotonically convergent algorithm Jl^l . In both cases the 
optimization was carried out with a standard penalty on the laser power, and the initial condition 
Pij(O) was set to thermal equilibrium at temperature T = 300 K. At this temperature the initial cis 
population is negligible. 

Below, the time-dependent population of the stable trans, is defined as the projection of p (t) 
onto the lowest 49 states localized around 0=0 and the cis population as the lowest 23 states 
localized around = it. The excited state population is the remainder. 
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FIG. 1: Time evolution under the optimized electric field with restriction on frequency and amplitude. 
Top panel: an electric field obtained by 64 iterations from random-initial phases ( . Middle panel: time- 
frequency resolved spectrum. Bottom panel: time evolution of populations. 

Figure [l] shows the isomerization dynamics under the optimized electric field restricted in fre- 
quency and amplitude. In the experiment |7|], the target was chosen as the ratio of the transmis- 
sion change AT at 400 nm to that at 460 nm, measured at 20 ps, assumed to be proportional to 
number of cis molecules created and number of trans molecules excited, respectively Q. The 
target is modeled computationally as the ratio of population of stable cis to the laser pulse area, 
/i ge / \E {t)\dt/h, at a target time 20 ps. Here the pulse area would be approximately proportional 

5 



to the number of trans molecules excited, if the nuclei were fixed and two electronic states par- 
ticipated in the process. Its use corrects somewhat for the presence of decoherence in this one 
dimensional model, relative to the multi-dimensional decoherence-free excitation experimentally. 
The electric field is optimized by 64 iterations from randomly initialized phases ; . The resul- 
tant electric field (Fig.[T]), in accord with experiment, is seen to have a peak at ~ 1.9 ps, and is 
considerably sharper than the laser 2 ps envelope for each frequency component, implying that the 
optimization has yielded a nearly equal set of phases 0,-. Figure[l]also shows small peaks around 
the main peak, giving an overall structure that is almost quantitatively the same as the experimen- 
tal result [7]. The bottom panel of Fig. \\\ shows the isomerization dynamics; the system is seen 
to be excited by the electric field at t ~ 2ps, with maximum P e (t) ~ 0.4. After photoexcitation, 
the system relaxes into stable trans and cis due to the system-bath coupling, giving a final cis 
probability of 0.17. The optimized electric field gives a target value of 0.53, while an unoptimized 
electric field gives a 0.7 times smaller value of 0.36, in agreement with experiment. 

Further studies shows that the small peaks in E(t) are of little relevance. That is, after 200 
iterations these small peaks almost disappear. Thus, the mechanism underlying the experimentally 
observed control is efficient photo-excitation under dissipative conditions, balancing the time scale 
for excitation and wave packet motion with the ongoing decoherence. 

Experimental results were also presented ED for the case of "no cis population", again using 
fields constrained in amplitude and frequency and using the transmission ratio as the target. The 
analogous computational result is shown in Fig. El in excellent agreement with experiment. Specif- 
ically, the pulse is now delocalized in time, rather complicated in form, and of duration longer than 
the dephasing time between the two electronic states. The maximum amplitude is considerably 
lower than that shown in Fig.[lJ Examination of the populations (bottom panel, Fig. O, shows that 
P e (t) reaches only 0.006, far lower than that in Fig.[T] Thus, despite the complexity of the pulse, 
its entire purpose is to ensure that there is no excitation of the initial trans species! However, the 
optimized electric field in our case gives a 0.05 target value, ~ 4 times smaller than the experi- 
ment, resulting from the fact that our model does not take into account competing processes such 
as excitations to other electronic states and the large number of nuclear vibrational modes. 

The robustness of our results were checked by performing various alternate computations. For 
example, using the target "ratio of created cis to depleted trans" was not useful since this quantity 
is always unity at 20ps. This is because in this model the molecule is, due to relaxation by 20ps, 
either stable trans or stable cis. Choosing a slower relaxation rate, so that excited population 



6 



10 



25400 



_ 25200 
§, 25000 
3 24800 
0.008 



£ 0.004 
& 0.002 
^ 0.000 



1 1 1 1 1 

1 1 1 1 1 


1 1 1 


1 


1 1 1 


PJt) 
- - p e (t) 


— / 

A i 


v ^1- - 



1 2 3 4 5 

f (ps) 



FIG. 2: Time evolution of the system under the optimized electric field that is restricted in frequency and 
amplitude. In this case, the target of the control is minimization of stable cis. Top panel: optimized electric 
field. Middle panel: time-frequency resolved spectrum. Bottom panel: time evolution of populations. 

remains at 20 ps, also did not improve the results. Here the target maximization showed an electric 
field moved to as early a time as possible to wait for the slow relaxation, or as late as possible 
in the case of minimization to avoid relaxation to trans and cis molecule. Neither agree with 
experiment. However, using of a simpler target, that is population of stable cis, gives qualitatively 
the same fields as shown in Figs [T] and El but dissimilar ratios of target improvement relative 
to the unmodulated pulse. Hence we are confident that the essential physics is contained in the 
simulation presented above. 

The experimental results optimize the target within the restricted frequency range and intensity 
described above. Ideally, however, adaptive feedback control desires the optimal result, which 
would require unrestricted laser equipment. To examine one such optimal solution we repeated 
the adaptive feedback studies with penalties on the power, but with no frequency restriction on the 
laser and with the population of stable cis used directly as the target. 

In the second of these target cases, the suppression of trans to cis isomerization, the optimal 
solution obtained was E(t) = 0, i.e. no excitation. By contrast, the fully optimal pulse for pro- 
ducing cis from trans is shown in Fig. where, for computational convenience, the target time 
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FIG. 3: Short-time evolution of the system under the fully optimized pulse. After 0.15ps, the electric 
field is essentially zero. Middle panel: time-frequency resolved spectrum of the pulse. Bottom panel: time 
evolution of populations. 

is 5 ps. The optimal electric field seen to consist of a pump pulse (from to 0.01 ps), and dump 
pulses (from 0.05 to 0.13 ps). The dump pulse has several frequency components, and is resonant 
with the deexcitation between electronic states around = 2.1 rad. After the pump pulse, P s (t) 
reaches almost 1, and after the dump pulse, ~ 0. 1 of population is transfered from P e (t) to P c is(t). 
Subsequently (not shown), the remaining excited component relaxes due to the system-bath cou- 
pling giving a final cis probability of 0.36, far higher than the cis probability obtained when this 
computation was repeated with the laser frequency restrictions above. 

Note that unlike the paradigmatic pump-dump coherent control scenario , or the pre- 

viously proposed eigenstate based cis-trans isomerization mechanism ifisfl the pump and dump 
steps are not coherently related since the system is decohered between pulses. Hence the optimal 
mechanism in this case is desirable excitation of an excited state population, delocalizing and de- 
cohering of the excited state wavepacket, and desirable de-excitation of the decohered excited state 
population. The delay between the pump and dump steps results from the spreading of density into 
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the cis region. 

It should be noted that the optimal pulse has almost zero amplitude after 0.15ps, since it is 
very difficult to try to make more cis after the dump process. That is, a second photo-excitation 
step from the stable trans to the excited state is not useful, since the excited component can not 
relax to stable cis within the few alloted ps. Also, due to the dephasing of the system, the excited 
component spreads rapidly into whole 2n range of <j>. Any attempt at a second photo-deexcitation 
from these excited component to stable cis would also be ineffective since it would be accompanied 
by compensatory excitation from the populated cis to the excited state. 

In summary, this work has (a) successfully exposed the simple underlying mechanisms associ- 
ated with the complex experimental results of an adaptive, condensed phase, feedback experiment, 
(b) demonstrated the role of frequency limitations in the experimental laser wavelengths and the 
concomitant emergence of an incoherent pump-dump scenario when these restrictions are lifted, 
with cis population as the target, and (c) demonstrated the utility of the simplest of models, one 
dimensional motion plus decoherence. The utility of such models derives from the fact that the 
decoherence is fast and that the measurement of isomer identity implicitly ignores all degrees of 
freedom but one (the angle <j>). Hence, one dimensional motion plus decoherence is formally the 
proper description. Note, however, that a appropriate representation of the decoherence is neces- 
sary to achieve the quality of results shown here. Since most liquid phase control experiments will 
be of a similar nature, simple control models of this type may be justified generally by the role 
and presence of decoherence. However, if decoherence is slower than the characteristic molecular 
dynamics then the problem is more complex and dynamics in many degrees of freedom must be 
explicitly controlled. 
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